CGPG-01/2-1 



Absence of Singularity in Loop Quantum Cosmology 



gr-qc/0102069 



Martin Bojc 

Center for Gravitational Physics and Geometry, The Pennsylvania State University, 
104 Davey Lab, University Park, PA 16802, USA 

It is shown that the cosmological singularity in isotropic minisuperspaces is naturally removed 
by quantum geometry. Already at the kinematical level, this is indicated by the fact that the inverse 
scale factor is represented by a bounded operator even though the classical quantity diverges at the 
initial singularity. The full demonstation comes from an analysis of quantum dynamics. Because of 
quantum geometry, the quantum evolution occurs in discrete time steps and does not break down 
when the volume becomes zero. Instead, space-time can be extended to a branch preceding the 
classical singularity independently of the matter coupled to the model. For large volume the correct 
semiclassical behavior is obtained. 
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On a macroscopic scale, the gravitational field is suc- 
cessfully described by general relativity, which is experi- 
mentally well tested in the weak field regime. However, 
this classical theory must break down in certain situa- 
tions where it predicts singularities, i.e. boundaries of 
space-time which can be reached by observers in finite 
proper time, but beyond which an extension of the space- 
time manifold is impossible [Q. An outstanding exam- 
ple is the big-bang singularity appearing in cosmologi- 
cal models. At this point curvature diverges whence the 
classical theory completely breaks down and has to be 
replaced by a quantum theory of gravity. However, up 
to now there is no complete quantum theory of gravity, 
and so the problem has been approached by first carry- 
ing out a symmetry reduction (by requiring isotropy and 
homogeneity) and then quantizing the resulting minisu- 
perspace models which have only a finite number of de- 
grees of freedom In the context of these models, as 
yet, there is no definitive resolution of the status of the 
initial singularity. Furthermore, generally the methods 
used in this analysis can easily miss some key features of 
the full theory. Indeed, while it has been speculated for 
a long time that quantum gravity may lead to a discrete 
structure of space and time which could cure classical sin- 
gularities, it has not been possible to embody this idea 
in standard quantum cosmological models. 

By now, there are promising candidates for a quan- 
tum theory of gravity. The results reported in this letter 
are obtained in the framework of quantum geometry ^] 
which does predict a discrete geometry because, e.g., the 
spectra of geometric operators as area and volume are 
discrete [^-Q]. Although temporal observables have not 
been included in the full theory, it is clear that the space- 
time structure is very different from that used in general 
relativity. But this difference can be important only at 
very short scales or in high curvature regimes like the 
one close to the classical singularity. This leads to the 
basic question raised here: What happens to the classical 
cosmological singularity in quantum geometry? 

The first step in our approach is the construction of 



isotropic states in full quantum geometry; we first quan- 
tize and then carry out a symmetry reduction. This, 
however, is not a straightforward problem because the 
discrete structure of space, represented by a graph (spin 
network) embedded in space, necessarily breaks any con- 
tinuous symmetry. But symmetric states can be defined 
as generalized states of quantum geometry which can 
be used for a reduction to minisuperspace models 
Note that this is not a standard symmetry reduction of 
the classical theory because symmetric states are inter- 
preted as generalized states in the full kinematical quan- 
tum theory. Only the Hamiltonian constraint has to be 
quantized and solved after the reduction. An immediate 
and striking consequence is that, in contrast to standard 
quantum cosmological models, spatial Riemannian ge- 
ometry is discrete leading to a discrete volume spectrum 
pO| . Furthermore, in contrast to standard quantum cos- 
mology, the same techniques as in the full theory can 
be used for the quantization of the reduced Hamiltonian 
constraint of the cosmological models [|l2|. This implies 
another difference, namely that the evolution equation is 
not a differential equation in time, but a difference equa- 
tion manifesting the discreteness of time p3|. 



Structure of isotropic models. According to 
states for isotropic models in the connection representa- 
tion are distributional states of the full kinematical quan- 
tum theory supported on isotropic connections of the 
form Al^ = cAjU!^ where A/ is an internal S'[/(2)-dreibein 
and oj^ are the left-invariant one-forms on the "transla- 
tional" part of the symmetry group acting on the space 
manifold S. The momenta are densitized triads of the 
form Ef — pAjXf with left-invariant densitized vector 
fields Xj fulfilling uj^{Xj) — Sj. Besides gauge freedom, 
there are only the two canonically conjugate variables 
{c,p} — K7/3 (k = SttG is the gravitational constant 
and 7 > the Barbero-Immirzi parameter) which have 
the physical meaning of extrinsic curvature and square 
of radius (a — ^/\p\ is the scale factor). The kinemati- 
cal Hilbert space Tikin = L'^{SU{2),Aiih) is the space of 
functions of isotropic connections which are square inte- 
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grable with respect to Haar measure. Orthonormal gauge 
invariant states are (see for details) 
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for j e |No together with (_i = (\/2sin|)"^ These 

states are eigenstates of the volume operator V with 
eigenvalues 110] 
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Later we will also use a different orthonormal basis of 
states adapted to the triad by introducing 
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where n represents the eigenvalues of p which determines 
the dreibein. In contrast to j, which is always positive 
and represents eigenvalues of the square of the scale fac- 
tor, n can also be negative. For this it is important that 
we have not only the character functions Xj i but also the 
additional functions Cj. This concludes the discussion of 
quantum states. 

The inverse scale factor. Classically, the metric of an 
isotropic spatial slice is given by qij = a^Sij = e^e^j 
where e\ is the co-triad. From this quantity we can build 
the expression 
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for the inverse scale factor, which we now quantize as a 
first application of the previously derived calculus. The 
co-triad is not a fundamental variable, but it can be quan- 
tized to 2i{'^lp)^^hi[h'^^ ^V] due to the classical identity 
ej, = 2(K7)-i{^a>^} 0- The expression dete in the 
denominator of to/j can be quantized to the volume op- 
erator which then can be absorbed into the commutators. 
Such a procedure has already been applied in in order 
to quantize matter Hamiltonians which become densely 
defined operators, and in the same way we arrive at the 
hounded operator 
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This operator is simultaneously diagonalizable with 
the volume operator and has the eigenvalues 
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where in the second step we have assumed that j — and 
hence Vj — is large. Thus, for large j, the leading term 
is the classical value V~^Sij, and the corrections (which 
are not necessarily isotropic) are of only the fourth or- 
der. We see that our quantization leading to a bounded 
operator does not spoil the classical limit. In fact, the 
a~^-behavior can be observed in a range which is much 
larger than expected from the large-j expansion. As Fig. 
|l| demonstrates, even for j — 1 are the eigenvalues very 
close to the classical expectation, and only the lowest 
three eigenvalues show large deviations. But this is al- 
ready deeply in the quantum regime, so such deviations 
are expected and lead to a finite behavior of the classi- 
cally diverging m//. Note that the volume operator has 
the eigenvalue zero (three-fold degenerate), but even in 
the corresponding eigenstates is the quantization of the 
inverse scale factor perfectly finite. This may be taken as 
a first indication for a removal of the classical singularity, 
although only at the kinematical level. 



FIG. 1. The classical expectation ^ (dashed line) and 
eigenvalues mii,j, j > of the inverse scale factor (x). Con- 
trary to the classical curve, the latter peak at j = 4 and 



decrease for j = and j ~ — \ {m 



is not shown). 



Discrete time evolution. Following the basic steps of 
the quantization in the full theory the Hamilto- 

nian constraints for cosmological models can be quan- 
tized with some adaptations to the symmetry For 
simplicity we write down here only the key term, the 
so-called Euclidean term H^^\ of the constraint opera- 
tor for spatially fiat isotropic models. However, all our 
qualitative results remain true for the full constraint and 
also for isotropic models with positive curvature. The 
constraint is of the form 



^^""^ = ^ E ^"''tTihlhjhj'hfhK[h],\V]) 
UK 

= sin^ I cos^ §(sin cos § - cos §V"sin |) 
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with action 

x(|n + 4) -2|n) + |n-4)). (6) 

In order to "unfreeze dynamics" and interpret solu- 
tions as "evolving states," as usual |jl^,|l6| we have to in- 
troduce an internal time which we choose as the dreibein 
coefficient p. Accordingly, we transform states |s) into 
an adapted representation by expanding |s) = J2n ^n\n) 
in eigenstates \n) oi p. This will allow us to find an in- 
terpretation of physical states as evolving histories. Fur- 
thermore, discrete geometry imples that eigenvalues of p 
are discrerte, whence time evolution is now discrete (see 
for details). Moreover, since we chose a geometrical 
quantity as time which can be negative and is zero for 
vanishing volume, we will be able to test the possibility 
of a quantum evolution through the classical singularity. 

To realize dynamics, we need to extend the model with 
matter degrees of freedom which can evolve with this in- 
ternal time. Matter can be incorporated by using co- 
efficients s„(0) depending on the matter field (f> in an 
appropriate fashion, the details of which is not impor- 
tant for what follows. The Hamiltonian constraint can 
then be written down using a matter Hamiltonian H,p 
(as in 10) which is diagonal in the gravitational degrees 
of freedom (and can also contain a cosmological term). 
The resulting quantum constraint equation can then be 
regarded as an evolution in discrete time: 

(^n+4|/2- "^n+4|/2-l)Sn+4('/')-2(Vj„|/2- V|„| /2_ i )s„ (0) 
+ (V\n-4\/2 - V|„_4|/2-l)s„-4(0) = jT'^^P S„(0) (7) 

{Vj are the eigenvalues (H) of the volume operator with 
VLi = 0) which is a difference equation for the coeffi- 
cients Sn(0) depending on the discrete label n (our dis- 
crete time). 

Fate of the singularity. Given initial data Sn{4>) for 
some negative n, we can use (Q) in order to determine 
later values for higher n. This, however, is possible 
only as long as the highest order coefficient T^„+4|/2 — 
V|„_i_4|/2-i is nonzero, which is the case if and only if 
n 7^ —4. So all coefficients for n < —4 are determined by 
the initial data. However, (^ does not determine so and 
instead leads to a consistency condition for the initial 
data. So the quantum evolution appears to break-down 
just at the classical singularity, i.e. at the zero eigenvalue 
of p. But this is not the case; in fact all s„ for n > 
are determined by (|^) from the initial data. This oc- 
curs because for n = we have: i) ^n|/2 ~ ^n|/2-i = 0, 
and ii) H^Sn{4') = 0; thus Sq completely drops out of 
the iterative evolution. E.g., s^ is determined solely by 
s_4 because the coefficient of s„ vanishes for n — Q. So 
we can evolve through the singularity and determine all 
s„ for n ^ 0. (The vanishing of H(j,sa{(f)) follows from 



the quantization of matter Hamiltonians ]14[ | similarly as 
described for the inverse scale factor.) 

Of course, in order to determine the complete state we 
also have to know sq, but a closer analysis reveals that 
So is fixed from the outset: The Hamiltonian constraint 
always has the eigenstate s„ = sqSuo with zero eigenvalue 
which is completely degenerate and not of physical inter- 
est. All evolving solutions are orthogonal to this state 
and have sq — which already fixes the coefficient sq left 
undetermined by using the evolution equation. We see 
that the complete state is determined by initial data for 
negative n, and so there is no singularity in isotropic loop 
quantum cosmology. The intuitive picture is as follows: 
Since for n < the volume eigenvalues V(|„|-i)/2 decrease 
with increasing n, there is a contracting branch for neg- 
ative n leading to a state of zero volume (in general, 
s±i 7^ and the volume vanishes for n = ±1 which cor- 
responds to j ~ 0) in which the universe bounces off lead- 
ing to the expanding branch for positive n which only can 
be seen in the classical theory and in standard quantum 
cosmology. This conclusion holds true for any kind of 
matter and cosmological constant, and is a purely quan- 
tum gravitational effect. In particular, we do not need to 
introduce matter violating energy conditions and thereby 
evade the singularity theorems. However, our result cru- 
cially depends on the factor ordering of the constraint 
which was chosen as one of the standard possibilities or- 
dering all triad components to the right. 

The semiclassical regime. We have seen that the clas- 
sical singularity is removed in loop quantum cosmol- 
ogy. But we need more for a viable cosmological model, 
namely we also need the correct behavior in the semi- 
classical regime. Classical behavior can only be present 
for large volume and small extrinsic curvature, i.e. if 
\n\ is large, c is small and the wave function does not 
vary strongly between successive times n (otherwise the 
state would have access to the Planck scale). In this 
regime we can interpolate between the discrete labels 
n and define a wave function ipi'^) ■= Sn(a)j n{a) := 
6a^ /jlp with a ranging over a continuous range (using 
a = y/\p\ ^ \/lh\/\n\/6 for large |n| as interpola- 
tion points). The difference operator A then becomes 
(As)„ := sn+i - Sn-i Ijlla'^dip/da + 0{ll/a^) 
leading to an approximate constraint operator i?*^^^ ~ 

-96(iA/2)2 . a/4 Q-f%{-ld/ii{a^)fa for large a. 

This is exactly what one obtains from the classical con- 
straint i?^^^ = —Qc^^/\p\ in standard quantum cosmol- 
ogy 101 quantizing 3c = —i^flpA/dp. In our frame- 
work, however, this is only an approximate equation valid 
for large scale factors. For this equation one can use 
WKB-techniques in order to derive the correct classical 
behavior. 

Going to smaller a one has to include more and more 
corrections in the expansion of the difference operators 
and also of the volume eigenvalues. By doing so one can 
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derive perturbative corrections for an effective Hamilto- 
nian including higher derivative terms. The closer we 
come to the classical singularity, the more corrections 
we have to include; and at the singularity we need to 
know all corrections which, as we know from our non- 
perturbative solution, have to add up to yield the dis- 
crete time behavior. So in these models higher order 
terms arise from the non-locality in discrete time of the 
fundamental theory. But even knowing all perturbative 
corrections, it would be very hard to see the correct be- 
havior without knowing the non-perturbative quantiza- 
tion. 

Quantum Euclidean space. In the simplest case, the 
Euclidean constraint for a spatially flat model without 
matter, it is possible to find an explicit solution to the 
constraint. The constraint equation is of order eight with 
one consistency condition as described above, so one ex- 
pects seven independent solutions. But we are interested 
only in solutions which have a classical regime in the pre- 
vious sense, i.e. no strong dependence on j for large j. 
Under this condition one can see that there is a unique 
(up to a constant factor) solution 



2j + l 
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in the connection representation. In standard quantum 
cosmology the constraint equation is (c) = with 

a solution ^/^[^(c) = (5(c) which is not unique. In or- 
der to compare the solutions we quantize a by axj — 
2t{jll,)-\V^+.-V^_i)Xj leading to ai^ cx E,(2j + l)Xj 
which in fact is the delta function on the configuration 
space SU{2). Therefore, we have a unique solution which 
incorporates the characterization of Euclidean space to 
have vanishing extrinsic curvature of its flat spatial slices. 

Conclusions. We have shown in this paper that 
canonical quantum gravity is well-suited to analyze the 
behavior close to the classical singularity. For this, it is 
important to use only techniques which are applicable in 
the full theory. This leads to a discrete structure of space 
and time which cannot be seen in standard quantum cos- 
mology. In our framework, the standard quantum osmo- 
logical description arises only as a limit for large volume 
where the discreteness is unimportant. For small vol- 
ume, quantum geometry leads to new effects which are 
responsible for the removal of the classical singularity. In 
contrast to earlier attempts this is not achieved by intro- 
ducing matter which violates energy conditions; it is a 
pure quantum gravity effect. It also does not avoid the 
zero volume state present in the classical singularity be- 
cause in general the wave function is not orthogonal to 
states with zero volume eigenvalue. Nevertheless there is 
no sign of a singularity because in quantum geometry it 
is possible to have vanishing volume but non-diverging 
inverse scale factor, which in isotropic models dictates 



all curvature blow-ups. Besides removing the singular- 
ity, the fact that an evolution through a state of zero 
volume is possible without problems could lead to topol- 
ogy change in quantum gravity. Technically, the removal 
of the singularity is achieved by using Thiemann's strat- 
egy |jll| of absorbing inverse powers of V into a Poisson 
bracket which also lead to densely defined matter Hamil- 



tonians 1 14 1 . So it is the same mechanism which regular- 
izes ultraviolet divergences in matter field theories and 
which removes the classical cosmological singularity. We 
have also seen that non-perturbative effects are solely 
responsible for this behavior and a purely perturbative 
analysis could not lead to these conclusions. 
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